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Using numerical simulations, we investigate the equilibrium dynamics of a single component fluid 
with Yukawa interaction potential. We show that, for a wide range of densities and temperatures, the 
dynamics of the system are in striking agreement with a simple model of generalized hydrodynamics. 
Since the Yukawa potential can describe the ion-ion interactions in a plasma, our results have 
significant applicability for both analyzing and interpreting the results of x-ray scattering data from 
high power lasers and fourth generation light sources. 
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I. INTRODUCTION 



Recently, using high power lasers and fourth generation 
x-ray sources, it has become possible to create and diag- 
nose extreme states of matter relevant to Inertial Con- 
finement fusion (ICF) and the cores of compact astro- 
physical objects in the laboratory (U-H)]. A particularly 
exciting development is that x-ray Thomson scattering 
experiments will soon be able to fully resolve time de- 
pendent ion dynamics in dense plasmas [H, H, Q • These 
ion dynamics are encoded in the wavevector and fre- 
quency dependent ion-ion structure factor (or simply dy- 
namical structure factor), Su(k,uj), which is the Fourier 
transform in space and time of the density autocorrela- 
tion function. For forthcoming experiments, an accurate 
model for the ion-ion structure factor is needed. 

In a previous work [§} , we found that the conventional 
hydrodynamic description (Navier-Stokes equations) re- 
produces Su(k,uj) well for k < k max , where k max X s ~ 
0.43 and X s is the electronic screening length. Despite 
the success of the conventional hydrodynamic descrip- 
tion at these large lcngthscales (small k), a model that 
works well at higher (momentum transfer) k is generally 
of greater applicability to the experiments. Fortunately, 
a well known framework - generalized hydrodynamics - 
already exists for extending the results of conventional 
hydrodynamics to these higher k values. In this paper, 
we compare one of the simplest models of generalized 
hydrodynamics to the results of state of the art numer- 
ical simulations for Sii(k,ui). We show that the model 
works remarkably well for all k values, i.e. the model 
describes both the conventional hydrodynamic limit at 
small k values and the large k behaviour (when the ions 
behave as a collection of free particles), along with the 
entire intermediate dynamics between these two regimes. 
Our results thus show that this simple model has sig- 



nificant applicability for analyzing and interpreting the 
results of forthcoming x-ray scattering experiments using 
fourth generation light sources. 

This paper is structured as follows. In Sec. [TlJ the 
Yukawa system - which represents interacting ions in a 
plasma - is introduced and details of our numerical sim- 
ulations of this system are given. In Sec. IIII1 the gen- 
eralized hydrodynamics framework is summarized, along 
with the Gaussian approximation for the memory func- 
tion that leads to a simple model for Su(k,Lu). This 
model is then shown to very accurately reproduce the 
results of our numerical simulations in Sec. IIV1 Also in 
this Section, we briefly discuss the applicability of our re- 
sults to x-ray scattering experiments (Sec. IIVE[) . before 
offering our conclusions in Sec. [V] 



II. NUMERICAL SIMULATIONS 

We consider a plasma consisting of one species of ions 
of charge Ze and mass m at temperature T and density 
n. Because the ions are much more massive than the 
electrons, on the time scale of the ion dynamics of in- 
terest here, electrons instantaneously screen the ion-ion 
Coulomb interactions and their degrees of freedom are 
not treated explicitly. We take the Yukawa potential, 



v(r) 



(Ze) 2 cxp(-r/A s ) 
47reor 



to represent the screened interaction between ions. The 
electronic screening length X s [3, [t| [lj| reduces to ei- 
ther the Debye-Huckcl law or the Thomas- Fermi distance 
in the limiting cases of classical and degenerate electron 
fluid respectively [l|. 

This single component system is known to be fully 
characterised by two dimensionless parameters only . 
These are: (i) the coupling strength 
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where a = (3/47rn) 1 / 3 is the average inter-particle dis- 
tance, and (ii) the screening parameter 



a 

aT 



In our MD simulations, we compute the dynamical 
structure factor, Su(k,w), of the Yukawa system for var- 
ious T values (1,5,10,50,120,175) at a = 0.1, 1.0 and 2.0, 
thereby spanning a range of thermodynamic conditions 
. In our simulations, the dynamics of N = 5000 par- 
ticles mutually interacting through the Yukawa poten- 
tial are resolved usin g th e Verlet algorithm in periodic 
boundary conditions [13]. In all cases, we include the 
Ewald summation in our force calculation - this is es- 
sential for small a values - using the particle-particle- 
particle-mesh (PPPM) method 13[. The rms error of 
our force calculation is 10 -5 . We find that obtaining 
accurate MD data for Su(k,ui) requires averaging the re- 
sults of a large number of simulations to improve statis- 
tics. This computational demand has made a thorough 
study such as ours impractical before now. For example, 
compared with the study of Hansen for the OCP sys- 
tem fbH ] - which, even after more than 30 years remains 
the primary source of MD data for quantitative studies 
of that system (l5| - we use 20 times as many particles, 
a smaller timestep by a factor of 2 — 10, and simula- 
tion times 200 — 40000 times as large. Ou r timestep 
St = O.OIcj" 1 , where u> p = \J (Z 2 e 2 n)/(eom) is the ion 
plasma frequency, ensures excellent energy conservation 
(AE/E m 10 -5 ). Moreover, we find that the long length 
of our simulations, 25 x 819.2 w" 1 for every T and a 
value, is of paramount importance: while it is possible 
to capture the essential features of Su(k,Lu) with simu- 
lations significantly shorter than this, producing a spec- 
trum that is of sufficient accuracy to draw conclusions 
about the validity of various models requires simulations 
of approximately this length (we note that our data for 
Sa(k,ui) changes negligibly by increasing the simulation 
time beyond 25 x 819.2 w" 1 ). In particular, these long 
simulation times are essential for computing accurately 
the decay time of collective modes at small k values (i.e. 
the width of the ion-acoustic peak in Su{k,u))). 

In a previous work Q, we presented MD results for 
Su(k,uj) of the Yukawa system at small k values; the 
MD data showed that the conventional hydrodynamic 
description works well in describing the dynamics pro- 
viding k < k max , where k max X s ~ 0.43. The new MD 
results presented here are for a significantly larger range 
of k values; in this paper we are interested in finding a 
model that reproduces the MD data for all k values. 



1 We have also performed some simulations at other a values; the 
model presented in Sec. IIII Al works very well for these other a 
values, but here we present results for a = 0.1, 1.0 and 2.0 only. 



III. MODEL 



A. Model for Su{k,u) 



In the hydrodynamic regime, the wavevector and fre- 
quency dependent ion-ion structure factor can be written 



s£(M = i {c s kfk 2 m 

S u (k) tt [w 2 - (c s k) 2 ] 2 + [wkhnf 



(1) 



where Su{k) is the static ion- ion structure factor. Equa- 
tion ([T|) is the result obtained from the linearised Navier 
Stokes equation [l2|, [l||. Here c s is the (isothermal) 
sound speed and rji is the kinematic viscosity. Equation 
([T]) clearly has considerable similarity to the expression 
that underlies the model we will consider in this article 



( w 2 )fcV(fc, w ) 



Sg(k,U)) _ 1 

S u (k) ~ n [u 2 - (w 2 ) - wfc V (k, w)] 2 + [wfc 2 0' (k, w)] 2 

(2) 

Equation ^ is a well known and exact representation 
of Sa(k,u>) that can be formally derived from micro- 
scopic theory 17[. The similarity to Eq. (p} is no co- 
incidence: Eq. ([2]) represents a generalized hydrodynam- 
ics in which both equilibrium properties and transport 
coefficients are replaced by suitably defined wavevector 
dependent quantities. In Eq. ©, (cj 2 ) = ^f- g 
defines a generalised isothermal sound speed c s (fc) = 

that, in the hydrodynamic 



limit of k — > 0, reduces to the conventional isothermal 

where xt is the 



sound speed c s (0) = c s = y^r^, 
isothermal compressibility of the system and Xt * na * °f 
an ideal gas. The quantities <j> (fc,w) and 4> (k,aj) are 
respectively the real and imaginary parts of the Laplace 
transform of the memory function <j>(k, t): in the analogy 
between Eqs. (fTJ) and ([2]), the memory function plays the 
role of a generalized viscosity. 

The model we present here amounts to using the Gaus- 
sian ansatz for the memory function, 



k 2 (j}(k,t) = k 2 (j>(k, 0) exp(-7rf 2 /4T^) 



| exp(-7rt 2 /4r j 



(3) 



where w 2 (fc) = (u A ) / (u 2 ) is given in terms of the fre- 
quency moments of Su(k,ui) 



uj n Su{k, u))duj . 



(4) 



Explicit expressions for (lo 2 ) and (ui 4 ) are given 

in the Appendix. Here Tfc, appearing in Eq. ([3|), is a 
wavevector dependent relaxation time. According to Eq. 
([3]), the real and imaginary parts of the Laplace transform 
of the memory function are given by, respectively (l4l[l8| . 



k 2 ^ {k,u) = [J > L {k)-((4)] Vc e 



7* 



(5) 



3 



and 

*?</>" {k,u>) = ^K(fc) - (ujDMt^/V^) , (6) 

where the Dawson function D{x) = 
cxp(-a; 2 ) f*exp(y 2 )dy 0. 

The quality of the Gaussian model has been previously 
identified for the Lennard-Jones fluid [ID, [2(| and by 
Hansen et al. in a pioneering study of the One Com- 
ponent Plasma (OCP) [l4j ; it has also been applied to 
experimental data for weakly coupled plasma produced 
by arc jets [2l[. However, because of the difficulty of 
conducting highly accurate numerical simulations at the 
time of the previous investigations, a detailed, conclusive 
comparison of the model in Eq. with the results of 
Molecular Dynamics (MD) simulations was not possible 
for those systems. Here, with the aid of modern com- 
puting facilities, we have conducted accurate, large scale 
MD simulations for Su(k,uj) across a wide range of ther- 
modynamic conditions. We find that the Gaussian model 
matches the MD data for the Yukawa system very well 
for all thermodynamic conditions we have examined in 
our simulations. 

B. Physical discussion of model for Su(k,u}) 



(O). Secondly, the viscosity is made to be non-local in 
space and time to model the viscoelastic effects in a real 
liquid. The memory function <p(r,t) that models these 
viscoelastic effects describes the delayed response of the 
longitudinal part of the stress tensor to a change in the 
rate of shear In Eq. ([3]), this response is mod- 

eled by a single relaxation time Tu- The requirement 
that the model reproduces the result obtained from the 
Navier-Stokes equations in the hydrodynamic limit gives 
a relation between the long wavelength behavior of this 
relaxation time and the kinematic viscosity rji [l8| , 

rn = mn lim>|(fc) - (iot)]r k /k 2 , (9) 

k— >0 

where r\i = (^ri + Q/mn, with r] and £ the shear and bulk 
viscosities respectively. 

The generalization included in Eq. ((5]) leads to the 
expression in Eq. © for the dynamical structure factor 
(see e.g. [H|). All that remains is to specify the memory 
function. As discussed in Sec. IIII Al here we choose a 
Gaussian memory function, as this is the simplest model 
that previous studies have suggested gives a good de- 
scription of the dynamics of classical fluids. We find that 
this choice yields a model of the dynamical structure fac- 
tor that matches the MD data for the Yukawa system 
remarkably well. 



The structure factor in the hydrodynamic regime, as 
given in Eq. ([I}, can be derived from the longitudinal 
component of the linearized Navicr Stokes equation, 

^- J(r, t) = --VP(r, t) + m V 2 J(r, t) , (7) 
at m 

where J(r,t) is the longitudinal current density and 
P(r,t) is the pressure. Similarly, Eq. @ can be de- 
rived from a generalized version of Eq. (JT)) (see [HI for 
more details), 

d T/ \ 1 „ f , > 5P(r,t) , . , . 

+ V 2 y f dsdr (f)(r - r ,t- s)J(r' ,s) , 

(8) 

where n(r, t) is the number density. This generalization 
is motivated in the following way. At small length scales, 
the validity of the conventional hydrodynamic descrip- 
tion can be expected to break down. Specifically, in the 
Navier Stokes description of Eq. ([7]), both the pressure 
term and viscosity term are local in space and time. The 
generalization in Eq. ([8]) includes the non-local behav- 
ior that is essential at small length scales in two ways. 
Firstly, it is assumed that a change in pressure at a po- 
sition r should not be determined completely by density 
fluctuations at the same position r but also by density 
fluctuations at neighbouring positions. This means that 
the pressure gradient due to a density gradient is non- 
local (hence the functional derivative appearing in Eq. 



IV. RESULTS AND ANALYSIS 

The Gaussian memory function model given in Eqs. 
d2]), (|5|) and (|6]) requires values for (w 2 ), u\{k) and for 
each k. Since all three of these parameters are in general 
unknown, we have fitted them to the MD spectrum of 
Sa{k,uj) using the least squares method. That is to say, 
for each k value for which we have computed Su{k,ui) 
with MD (these are the k values compatible with the 
periodic boundary conditions in our simulations), we fit 
the model to the MD spectrum of Su(k, to). When this is 
done, the model reproduces the MD data very accurately 
for all r and a values; in Sec. IIV Al we show that this is 
the case for small, intermediate and large k values (see 
also [H|)- 

The three parameter fit is the correct way to compare 
the Gaussian memory function model to the MD spec- 
trum of Su{k,uj). This is true despite the fact that two 
of the parameters, (cjy and wf^fc), can in principle be 
obtained by computing Sn (k) (or equivalently the radial 
distribution function g(r) [12| |) with MD and using the 
formulae given in the Appendix „ When obtained from 
MD in this way, these two parameters are subject to nu- 
merical incertainty. Therefore, one would expect that 
constraining (w 2 .) and w 2 (fe) - and therefore fitting the 
model to the MD spectrum using only a single parame- 
ter Tfc [13, [H, HU - would result in poorer fits and larger 
errors. In Fig [TJ we show that in general this is indeed 
the case. 

The validity of the three parameter fit can be confirmed 
by comparing the fitted values of the two parameters (u> 2 ) 
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FIG. 1: (color online) Comparison between the Gaussian 
model when only the parameter is fitted to the MD spec- 
trum (dashed line), and when all three parameters are fitted 
(solid line) for four separate cases. The MD results are given 
by the dots. 



and oj\(k) to their values when instead computed with 
MD as described above. As shown in Figs. [5] and [31 the 
parameters and oj\{k) obtained from the fit to the 
MD spectrum of Su{k,ui) agree very well (within 10%) 
with those computed from the MD g(r) and Su(k). This 
is only the case because the model works so well. For 
example, as shown in Fig. [3j if an exponential rather 
than Gaussian memory function is used (this is known 
as the viscoelastic model and is discussed in Sec. IIV CI) , 
the numerical values obtained for wf^fc) by fitting the 
model with three parameters do not agree well with those 
computed from the MD g(r) and Sa{k). In the remainder 
of the paper, we present only the results for the Gaussian 
memory function model with three fitting parameters; 
the one parameter fits are irrelevant as their comparison 
with the MD data for Sa(k,u>) is not indicative of the 
quality of the model. 



Comparison between model and MD 
simulations 



We find that in general the Gaussian memory func- 
tion model reproduces the MD data very well for all of 
the r (1,5,10,50,120,175) and a (0.1,1 and 2) values we 
have considered, at all k values (our simulations are for 
ka = 0.23 — 6.19). Extended figures of our complete MD 
results are available as supplementary material [22J ; here, 
in Figs. 0] -El we show only a selection of these complete 
results at small, intermediate, and large k respectively. 

At small k values (Fig. 2]), for all a and T, the MD data 
shows a clear ion-acoustic (or Brillouin) peak that repre- 
sents a damped sound wave in the plasma. In this regime, 
the model extends the conventional hydrodynamic de- 
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FIG. 2: (color online) Comparison between (u;^) as computed 
from MD using the formulae in the Appendix (dashed line, 
with 10% error band), and the values obtained from the three 
parameter fit of the Gaussian memory function model (trian- 
gles) and the viscoelastic model (squares) for three different 
plasma conditions, (a) T — 120, a = 0.1, (b) T = 50, a = 1, 
(c) T = 175, a = 1. 



scription to finite k values. Specifically, the generalised 
sound speed along with the imaginary part of 4>(k, lo) cor- 
rect for the fact that the position of the peak does not 
vary linearly with k as in the hydrodynamic description 
[H, and the real part of 0(fc,w) corrects for the width. 

At intermediate k values (Fig. [SJ, the model gives a 
surprisingly accurate account of both the width and posi- 
tion of the ion acoustic peak. This is particularly true for 
r < 50. For higher T values, the MD data does in some 
cases show additional structure which the model cannot 
recreate. In particular, for a = 0.1 and 1, a two peak 
structure is visible for ka = 2.32 and a three peak struc- 
ture for ka = 3.09 (e.g. Fig. [51 top left). The small peak 
just below uip for ka — 3.09 is of particular interest - it 
does not appear to have been seen or commented upon in 
previous MD calculations. We note that this peak is dis- 
tinct from the higher harmonic peaks reported in [23j . In 
fact, at k = 0.1 only, we do see signs of a second harmonic 
peak, at a frequency close to 2w p . We have neglected this 
harmonic peak in our analysis, since we find it to be more 
than than 3 orders of magnitude smaller than the main 
features in the spectrum of Su(k,uj), in good agreement 
with [23[. On the other hand, the peak shown in Fig. [5J 
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FIG. 3: (color online) Comparison between (wf,(fc)) as com- 
puted from MD using the formulae in the Appendix (dashed 
line, with 10% error band), and the values obtained from the 
three parameter fit of the Gaussian memory function model 
(triangles) and the viscoelastic model (squares) for three dif- 
ferent plasma conditions, (a) F = 120, a — 0.1, (b) V = 50, 
a = 1, (c) T = 175, a = 1. 
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FIG. 4: (color online) Comparison between the MD data for 
Su(k,Lu) (dots) and the Gaussian memory function model 
with three fitting parameters (solid line) for small ka values. 



(top left) is of the same order of magnitude as the main 
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FIG. 5: (color online) Comparison between the MD data for 
Sa(k,oj) (dots) and the Gaussian memory function model 
with three fitting parameters (solid line) for intermediate ka 
values. 
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FIG. 6: (color online) Comparison between the MD data for 
Sa(k,oj) (dots) and the Gaussian memory function model 
with three fitting parameters (solid line) for large ka values. 



features of Su(k,u>). We believe that this peak is due 
to microscopic 'caging' effects (e.g. [TU, Qj})- That is, 
at these lengthscales, the relatively high frequency oscil- 
lations of individual particles in the cages produced by 
their neighbors are imprinted on Su(k,io). We note that 
although the model does not fully capture the additional 
structure in the MD data for these conditions, on average 
it does give a good account of the overall shape of the 
spectrum. 

At large k values (Fig. [5]), Su(k,uj) reduces to a single 
peak at lo = 0. In this regime, the model reproduces 
the MD data very accurately in all cases. As k increases, 
Sa(k, ui) should tend to its ideal gas limit S^(fc, w), which 



6 



is independent of a [IH, [l4| , 

^> = fef-(-^)- (10) 

As shown in Fig. [7J at constant a, as T increases Su(k, u) 
converges more slowly towards S^k,^). Indeed, at the 
highest k value we have considered in our MD simulations 
(ka = 6.19), the MD result only compares well to its 
ideal gas limit for T < 10 (see Fig. [7J. We note that 
the discrepancy between Su(k,Lu) and its ideal gas limit 
can more readily be seen by looking at the MD data for 
the static structure factor Su{k) ; the ideal gas limit will 
only be approximated at k values for which Su(k) ~ 1 
(since Sg(ft) = 1). 
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FIG. 7: (color online) Comparison between the MD data for 
Su(k,uj) for a = 1 and ka = 6.19 (dots) and the ideal gas 
limit given by Eq. (|10|) (solid line). Also shown on each panel 
is the value of the static structure factor Sa(k) at ka = 6.19. 

In any case, as shown in Fig. [51 the Gaussian model 
compares very well to the MD data at our highest k value 
of ka = 6.19, regardless of whether or not this k value is 
sufficiently large for Sn(k,u>) to be close to its ideal gas 
limit. 



B. Hydrodynamic limit 

In previous investigations (e.g. (3), Eq. ^ was 
used to infer the kinematic viscosity from the long wave- 
length behavior of the relaxation time appearing in the 
memory function. For the Yukawa system, in principle 
this could be used to determine the shear viscosity (the 
bulk viscosity is in general negligible in comparison with 
the shear viscosity for the Yukawa system [Hj]). How- 
ever, due to the inaccuracy inherent in measuring the 
width of the (very narrow) ion acoustic peak obtained 
from the MD simulations at small k values, we find that 



this method is of little practical use compared to other 
approaches to determining the viscosity. These alterna- 
tive approaches include utilizing the Green-Kubo relation 
for the shear stress autocorrelation function (25|, non- 
equilibrium molecular dynamics methods [26| . and com- 
putation of the transverse current autocorrelation func- 
tion [13. 

Along with the generalized viscosity, as discussed in 
Sec. IIII Al in the hydrodynami c limit k —> the gen- 
eralized sound speed c s (k) = \/ (wf) / k 2 reduces to the 
conventional (isothermal) sound speed c s . The small k 
behaviour of the generalized viscosity and sound speed 
thus ensure that using the Gaussian ansatz for the mem- 
ory function in Eq. ([2]) gives a result that is compati- 
ble with the result obtained from the linearised Navicr 
Stokes equations fl2l ] when thermal fluctuations are ne- 
glected. To be clear, Eq. is an entirely general (i.e. 
exact) representation of Su(k,u)), The effective neglect 
of thermal fluctuations is made by assuming the ansatz in 
Eq. ([3]). That is to say, in the case of the Gaussian ansatz 
it is instructive to think of the memory function as a sort 
of generalized viscosity. There is no term in the memory 
function that represents the effects of temperature fluc- 
tuations i.e. a generalized (or indeed non-generalized) 
thermal conductivity. 

It is straightforward to modify Eq. ([3]) so that the 
result from the Navier Stokes equations including tem- 
perature fluctuations is recovered in the hydrodynamic 
limit (see e.g. [l(| tH)]). The simplest extension in- 
volves maintaining a generalized sound speed and vis- 
cosity, and adding the (non-generalized) thermal con- 
ductivity contribution obtained from conventional hy- 
drodynamics (the Navier-Stokes equations) as an addi- 
tional term in the memory function. In a more involved 
scheme, this additional contribution can also be general- 
ized dint 

For the Yukawa system with the T and a values we 
have considered here, including in the memory function 
the effects of thermal fluctuations is unnecessary. This is 
because the ratio of specific heats, 7, is very close to 1, 
as indicated by the absence of a Rayleigh peak at U) = 
for small k in the MD data (Fi g. EQ ), as well as previous 
equation of state calculations [28|]. The only cases in 
which this peak - which represents a diffusive thermal 
mode - is not negligible is for the more weakly coupled 
(r < 10) systems at a = 2 (see Fig. [U bottom left). As 
expected, the model does not capture this peak in the 
MD data. 

The fact that 7 w 1 for the Yukawa system with the 
r and a values considered here is certainly a reason why 
the Gaussian memory function works so well. Indeed, 
the ansatz in Eq. ^ would not be expected to work 
as well when the ratio of specific heats 7 is noticeably 
different from unity [l~7| ; this includes the Yukawa system 

for r< 1. 
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C. Comparison with viscoelastic model 

Given the excellent agreement between the MD data 
and the Gaussian memory function model, we have not 
found it necessary to undertake an exhaustive compari- 
son with the numerous other forms of memory function 
proposed in the literature [l6j . However, here we briefly 
comment on another widely studied and used ansatz for 
the memory function 



k 2 cf>(k,t) 



k 2 <f>(k,0)exp(-t/T%) 
[ul{k)~{ul)] C M-t/^) 



(11) 



When combined with Eq. @, Eq. (fTTj) - which repre- 
sents the simplest assumption that can be made about 
the time dependence of the memory function - is known 
as the viscoelastic model [l3] ■ 

As indicated in Fig. [5] and discussed in detail else- 
where [13, EM , the viscoelastic model cannot capture 
the shape of Sa(k,ui) across a large range of k values. 
While the model works well at small k (indeed, for the 
viscoelastic model the results of isothermal hydrodynam- 
ics are again recovered, with a relation between the re- 
laxation time and the kinematic viscosity similar to 
Eq. (O), the model tends to predict rather more struc- 
ture in Sii{k,uS) than is evident in the MD data (Fig. 
|8]). Clearly then the Gaussian memory function is vastly 
superior to the exponential one. 
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FIG. 8: (color online) A sample of our MD results for Su(k, w) 
at r = 10, a = 1 (dots) contrasting the results of the model 
in Eq. ([2]) for exponential (dashed line) and Gaussian (solid 
line) memory functions. 
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FIG. 9: (color online) The relaxation time tu as determined 
from the fit of the Gaussian model to the MD spectrum of 
Su(k, uj) for a = 2 and a range of F values. 



of Sa(k,uj) for a = 2. As shown in Fig. [9j we find 
that as k increases, Tk decreases. This agrees qualita- 
tively with e.g. the behavior of the relaxation time de- 
termined for the Lennard- Jones fluid in previous investi- 
gations [ill, [29[ . One certainly expects that at decreas- 
ing wavevectors, the relaxation time should increase: as 
k — > 0, the memory function should decay fast enough to 
guarantee the validity of the Markovian approximation, 
which itself is related to the fulfillment of the conserva- 
tion laws (l7j . 

In our investigation, we find that at the very smallest 
k values accessible to our simulations (i.e. below ka = 
0.64, which is the minimum k value shown in Fig. [5]), 
the numerical value of Tfe is difficult to extract from the 
MD spectrum reliably, and therefore it is not possible 
to examine the exact k — > behavior of the relaxation 
time. That is to say, the fitted value of Tk at these small k 
values does not connect smoothly to the values at higher 
k values; this is because the spectrum Su(k,uj) consists 
of a very sharp peak, for which it is difficult to accurately 
determine the parameters in the Gaussian model (see also 
Sec. ITVBl . 

Physically, the relaxation time Tk controls the specific 
collective behavior of the system: for times t <C Tk the 
system responds 'elastically' (i.e. like a 'frozen' solid-like 
system), wheras for times t ^ Tk the viscous mecha- 
nisms set in and reveal the inherent dynamic disorder 
[J3]- Therefore, the decrease in Tk as k increases corre- 
sponds physically to the fact that at increasingly short 
lcngthscales, viscous behavior is observed at increasingly 
short timescales. 



E. Applicability to x-ray scattering experiments 



D. Discussion of the relaxation time Tk 

Fig. [5] shows the relaxation time Tk as determined 
from the fit of the Gaussian model to the MD spectrum 



In a previous work [§J, it was shown that the conven- 
tional hydrodynamic description (i.e. Eq. ([!])) is valid 



providing k < k 



where k max X s ~ 0.43. This means 



that experiments designed to measure Su(k,oj) [l[ at k 
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values below k max can in principle be used to determine 
transport (e.g. viscosity) and thermodynamic properties 
(e.g. compressibility) of dense plasmas. 

At k values larger than k max , our results show that 
the Gaussian memory function model extends the con- 
ventional hydrodynamic description very satisfactorily. 
Thus experiments for k > k max measure the generalized 
quantities appearing in the memory function model of 
Eq. ©. 

Present x-ray scattering experiments are also con- 
cerned with diagnosing the density and temperature of 
dense plasmas For this task theoretical models for 
how Su(k,uj) depends on density and temperature are 
required. In the Yukawa system, the density and temper- 
ature are encoded in T and a. Thus here we briefly look 
qualitatively at how Su(k, uj) changes with T and a: this 
gives an indication of how the experimental scattering 
cross section should vary with density and temperature. 
We restrict ourselves to the region of k values for which 
S'ii(fc,w) shows a clear ion-acoustic peak, since then its 
description reduces to the position, width and height of 
this peak. 

Fig. [TOlshows how the position, width and height of the 
ion-acoustic peak as extracted from our MD simulations 
vary with reduced wavenumber ka for a number of F and 
a values. As shown in the top left panel of Fig. [TU1 
the dependence on ka of the ion-acoustic peak position 
is almost identical for a large range of T values (i.e. T = 
50 — 175). The peak width and height do show more 
discernible differences for these T values. At smaller T 
values (r = 1,5 and 10), the differences in the position, 
width and height of the peak are greater. 

At constant T (right panels of Fig. [TU)) . the peak posi- 
tion is rather different for a = 0.1,1 and 2. In this case, 
the width and height are more similar, particularly for 
a = 0.1 and a = 1. 

We expect that a given experiment will be able to de- 
termine peak position, width and height at a specific 
wavenumber (determined by the scattering angle and x- 
ray wavelength [l[). The extraction of T and a values 
could then be done by using these experimental results 
in conjunction with a set of three plots as shown in Fig. 

EH 

Of course, our discussion in this section assumes that 
a real physical plasma at a certain density, temperature 
and (average) ionization state can be described by the 
Yukawa system. While in principle this mapping could 
be attempted for any given values of these plasma param- 
eters, our main interest at present concerns the dense (ap- 
proximately solid density), liquid-like plasmas at temper- 
atures of « lOeV that can be created in high power laser 
experiments [l[. Recently, a method for mapping the 
physical parameters of these states to the Yukawa model 
(i.e. determination of T and a) has been suggested [3fj| . 
Therefore, we expect that the results we have obtained 
for the Yukawa system are certainly relevant for future 
experiments that will measure ion dynamics of these ex- 
treme states of matter. 
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scale) and height (log scale) of ion-acoustic peak with reduced 
wavenumber ka. The left panels are for a — 1 and a range of 
F values, and the right panels are for T = 10 and a range of 
a values. 



V. CONCLUDING COMMENTS 

The Gaussian memory function model is an extremely 
good representation of the dynamical structure factor 
Su(k,uj) of the Yukawa system for a wide range of ther- 
modynamic conditions. The model very accurately re- 
produces the spectrum of Sn{k,u>) from MD in terms 
of just 3 parameters and, as such, it is a useful way of 
accurately condensing or representing such data. This 
conclusion was only possible because of the highly accu- 
rate MD data presented in this paper. The model can be 
used by fitting either a single parameter or three param- 
eters to the spectrum of Sn{k,w) at a particular k value; 
in the latter case, the small numerical inaccuracies that 
arise in the MD simulations can be accounted for. 

Why exactly this form of memory function should work 
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so well is an interesting question that certainly merits 
further investigation. Other memory function models, 
such as the viscoclastic model (an exponential memory 
function) do not compare well to the MD data for a wide 
range of k values. It is possible that the reason a faster 
decaying (compared to exponential) Gaussian works well 
is related to the chaotic nature of classical systems - this 
is reflected in the relatively short 'memory' of the system. 

Since the Yukawa system can describe ion-ion inter- 
actions in a plasma, our results are applicable to future 
x-ray scattering experiments that will attempt to mea- 
sure ion dynamics in dense plasmas 0- In particular, our 
MD results for the position, width and height of the ion- 
acoustic peak could be used to infer the thermodynamic 
conditions of dense plasmas. 
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are given in terms of the frequency moments of Sa(k,ui), 
defined as 



(u n ) = / w n S u (k,u)du. 



(A.3) 



The zeroth moment of Su(k,cj) gives the static structure 
factor Sii(k) 



The second moment is 



U7, 



Sa(k) . 



3r 



(A.4) 



(A.5) 



where q = ka is the reduced wavevector (a = 
(3 /(47m)) 1 / 3 is the Wigner-Seitz radius) and uj p = 
y/ (Z 2 e 2 n) I (com) is the (ion) plasma frequency. The 
fourth moment is (see 17J, Eq. (1.137)) 



1 

3T 



r 



+ q 2 n% - q 2 M(qf, af) 



(A.6) 



Here r = r/a, the Einstein frequency VLe is given by 

,2 roc 



n 2 



r exp(— ar)g(r)dr , (A. 7) 



Appendix: Frequency moments of Su{k,ui) 
The wavevector dependent quantities, 



and 



ujl(k) 



M "' ~ (lo 2 ) 



(A.l) 



(A.2) 



and 

M(x,y) 



1 



g(r) exp(-y) 



I y +2/+1 ) x 



sin x 3 cos x 3 sin x 



y aw 
3a; 



dr . 
(A. 



Eqs. IA.6I - lAJjl give an exact expression for the fourth 
moment for the Yukawa one component plasma. 
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